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ABSTRACT 

We explore the cosmological implications of 7 deep survey fields observed at the Berkeley-Illinois- 
Maryland Association (BIMA) Array with 30 GHz receivers. These observations probe the Cosmic Mi- 
crowave Background anisotropy on scales corresponding to I ~ 5500, and an earlier analysis of these data 
detected no galaxy clusters. We use numerical cluster simulations and mock observations to characterize 
the cluster detection efficiency for each of the BIMA fields. With these detection efficiencies we derive 
constraints on the cosmological parmaters Q m and <7g > ruling out those models which overproduce galaxy 
clusters. Using only these seven BIMA fields, we calculate a 2a upper limit of a 8 < I.00^°- 43nM ~ ' 22 
for flat models with 0.1 < Q.m < 1- When the power spectrum of density fluctuations is COBE normal- 
ized, we find VIm < 0.63 at 95% confidence level for flat models. This constraint includes our estimate of 
the large uncertainties in the SZE flux-virial mass relationship as well as published uncertainties in the 
Hubble parameter, the COBE power spectrum normalization and the primordial power spectrum index. 
In addition, we account for the effects of sample variance. Thus, we conclude that the non-detections 
are to be expected in a low Qm universe, given the sensitivity and solid angle of the deepest SZE survey 
to date. 

Subject headings: cosmology: observation - galaxies: clusters - cosmic microwave background 



scale CMB anisotropy (Holzapfel et al. 2000a). Contin- 
ued observations on these and additional fields has led to 
a detection of anisotropy with flat band power at the level 
Qfiat — 6.li 2 'g/i-K", where uncertainties describe 68% con- 
fidence regions (Dawson et al. 2001, see also Dawson et al. 
2002). In addition, no galaxy cluster detections were re- 
ported, which is interesting because of the sensitivity of 
SZE surveys to high redshift clusters. In the currently fa- 
vored cosmology, these seven fields sweep out as much co- 
moving volume to z ~ 2 as a 130 deg 2 local survey sweeps 
out to z ~ 0.1. 

In principle, cluster surveys constrain cosmological pa- 
rameters through the effects those parameters have on the 
volume surveyed per solid angle and the evolution of clus- 
ter abundance (e.g. Haiman et al. 2001; Holder et al. 2001; 
Mohr 2001). In this analysis we focus on the lack of clus- 
ters in these deep '97 and '98 BIMA fields. A less sensitive, 
earlier survey with the Australian Telescope Compact Ar- 
ray also provided upper limits on the CMB anisotropy on 
similarly small scales (Subrahmanyan et al. 2000). An 
analysis of the cosmological implications of that survey, 
which relied on the assumption that the dominant source 
of anisotropy on arcminute scales is the thermal SZE from 
clusters, concluded that a Iow-JIm universe is preferred 
(Majumdar & Subrahmanyan 2000). 

An analysis of this BIMA survey allows us to rule out 
cosmological models which produce too many clusters. 
The expected number of detected clusters in a survey is 

(N) = j dnj dV J™ dM f(M) £L (M, z ) , (1) 

where (dn/dM)dM is the comoving number density of ob- 
jects of mass between M and M + dM, f(M) is the detec- 
tion efficiency (0 < f(M) < 1), which takes into account 
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1. INTRODUCTION 

Cosmic microwave background (CMB) images provide 
information about the presence of galaxy clusters over 
a wide range of redshift, because galaxy clusters leave 
signatures in the CMB through the so-called Sunyaev- 
Zel'dovich effect (SZE; Sunyaev & Zeldovich 1970; Sun- 
yaev & Zel'dovich 1972). The SZE is produced by inverse 
Compton interactions between the CMB photons and the 
hot electrons (fcsT e ~ IkeV, see §2.1) in the intracluster 
medium. The intracluster medium is shocked and heated 
in the process of cluster formation. The net result in the 
SZE is a transfer of energy from the electron population 
to the CMB photons, leading to a distortion in the CMB 
spectrum. The thermal SZE signature of a cluster in 1 cm 
(~ 30 GHz) observations is a reduction in the CMB bright- 
ness. On angular scales of ~ 5 and smaller the thermal 
SZE contribution to the CMB anisotropy can dominate 
that of the primary CMB anisotropy (Holder & Carlstrom 
1999; Holzapfel et al. 2000a; Hu & Dodelson 2001; Springel 
et al. 2001). 

Here we discuss the cosmological implications of 7 deep 
survey fields observed at the Berkeley-IUinois-Maryland 
Association (BIMA) array with the 30 GHz receivers 
(Carlstrom et al. 1996; Holzapfel et al. 2000a). The ob- 
servations were made during 1997 and 1998, in a com- 
pact configuration at 28.5 GHz, providing a Gaussian pri- 
mary beam with FWHM~ 6'6. This survey, originally 
planned to probe the CMB anisotropy on arcminute scales 
I ~ 5500 (corresponding to angular size ~ 2'), provides 
relatively large sky coverage (~ 250 arcminute 2 ) for an 
arcminute-scale anisotropy experiment. These observa- 
tions have been used to place an upper limit on small 
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the fact that only clusters with masses above some mass 
limits can be detected. The first integral is over the solid 
angle of the survey, the second integral is over the redshift, 
and the third integral sums over the portion of the cluster 
mass function to which our survey is sensitive. In addi- 
tion to its dependence on mass, the detection efficiency 
f(M) depends on the sky position (a, S), the redshift z 
and cosmological parameters such as H , Am and 

Cosmological constraints come through comparison of 
the observed and the expected number of clusters. In the 
BIMA survey no clusters were detected; assuming Pois- 
son statistics, the probability P of observing clusters is 
P(0\(N)) = e _ W. Using this approach we explore the pa- 
rameter space spanned by Qm and <j%, checking the con- 
sistency of these cosmological models. We also consider 
the uncertainties in the cluster detection efficiency f(M), 
Hubble parameter, COBE power spectrum normalization, 
power spectrum index and the effects of sample variance. 

This paper is organized as follows: in §2 we present 
a general description of interferometric SZE observations, 
describe the BIMA survey data, and determine the sur- 
vey detection efficiency via mock observations of hydro- 
dynamic cluster simulations. We discuss the ingredients 
needed for the calculations of cluster abundance in §3. 
The cosmological constraints and a discussion on effects 
of systematic errors appear in §4. We present a summary 
and discussion on prospects of future interferometry SZE 
observations in §5. We discuss the Ho scaling of interfer- 
ometric SZE observations in the Appendix. 

2. SURVEY SENSITIVITY TO CLUSTERS 

The expected number of detected clusters in each cos- 
mology depends sensitively on the detection efficiency 
function f(M), characterized by a limiting detectable mass 
Mu m at which f(Mu m ) — 0.5. We determine the detec- 
tion efficiency by conducting mock observations of hydro- 
dynamical cluster simulations using the exact character- 
istics of the real BIMA observations. In this section, we 
first discuss the interferometric observations of the cluster 
SZE, then describe the BIMA data, and provide a detailed 
description of the determination of the detection efficiency. 

2.1. Interferometric SZE Observations of Clusters 

Interferometric observations of cluster SZE have greatly 
advanced in the past decade (Carlstrom et al. 2000). In- 
terferometers can observe a wide range of angular scales si- 
multaneously, reduce possible confusion due to radio point 
sources toward galaxy clusters, and generate high-fidelity 
SZE images. 

Interferometers measure the visibility V\(u,v), whose 
Fourier transform pair, the specific intensity I\(9 x ,9 y ), is 
directly proportional to the Compton y(6 Xl 6 y ) distribu- 
tion of the cluster. Along any line of sight, 

y = a T [ dlnJ^-\, (2) 
J m e c 2 

where ctt is the Thomson scattering cross-section, n e is the 
electron number density, m e is the electron rest mass, c is 
the speed of light, ks is the Boltzmann constant and T e 
is the electron temperature. The SZE temperature decre- 
ment is then AT = g(x)yTcMB, where x — hv /ksTcMB 
is the dimensionless frequency, Tqmb = 2.728 K (Fixsen 



et al. 1996) is the mean cosmic microwave background 
temperature, and, within the Kompaneets approximation, 
g(x) = x(e x + l)/(e x — 1) — 4 is the amplitude of the spec- 
tral distortion. In the Rayleigh-Jeans regime (x « I) 
g(x) s=s —2, and at the central frequency of these BIMA 
observations, v = 28.5 GHz, g(x) = —1.9583. We con- 
vert this decrement distribution AT(9 x ,0 y ) into specific 
intensity I\(9 x ,9 y ) using 

h(e , M = 2 ^^M. (3) 

where A is the wavelength of observation. The visibility 
V\(u, v) then is the Fourier transform of the intensity dis- 
tribution or image I\(9 x ,9 y ) multiplied by the telescope 
primary beam B(9 x ,9 y ), 

V x (u,v) = J J d9 x d9 y / A B e - 2 «H.W„) 

= I~x ® B, (4) 

where I\ and B are the Fourier transforms of I\ and B, 
respectively, and (g> denotes convolution. For this analysis, 
we take the BIMA primary beam to be a Gaussian with 
FWHM= 6'. 33; therefore B is also a Gaussian. Because 
Va(0) = J d 2 9I\B, the overall normalization of the mea- 
sured visibilities in the case that the cluster lies well within 
the BIMA primary beam (i.e. B ~ I ) and observation fre- 
quencies are in the Rayleigh-Jeans limit is directly related 
to the cluster flux within the virial region: 

T/ / n N b~i a -4 a T k B T C mb /rp x fiCM M 2 oo , , 

v x (0) — s 200 = ^_^^<T e > — — — (5) 

where (T e ) is the electron density weighted mean tempera- 
ture of the cluster, cU is the angular diameter distance, /i e 
is the mean molecular weight per electron, m p is the rest 
mass of proton, ficM is the ratio of intracluster gas mass 
to binding mass, and M 2 oo = 47r/3r 200 x 200p c is a mea- 
sure of total cluster binding mass, within which the mean 
overdensity is 200 times of the critical density p c . This ex- 
pression also defines the virial radius r 20 o • Note that Eqn 5 
does not account for any SZE flux contribution from out- 
side the cluster virial region. Also notice that, the only 
cosmology-dependence of the flux comes from dA- 

To examine the sensitivity of cluster visibilities to clus- 
ter size and morphology, we adopt a cluster tempera- 
ture decrement model. We use the spherical /3 model 
(Cavaliere & Fusco-Femiano 1978), AT (6) = AT (1 + 
(9 / 9 c ) 2 y i ~ 3/3 ^ 2 , where AT is the temperature decrement 
toward the cluster center, and 9 C — r c /dA is the angular 
core radius. Comparison of this model with our simulated 
cluster AT maps indicates that it provides a reasonable 
description of clusters with (/?) ~ 1.1. For (3 = 4/3 in the 
case where the primary beam is larger than the cluster (i.e. 
B ~ 1), the visibility and the SZE flux within the virial 
radius is related by: 

V x (q) = ^e" 2 ^ (6) 

where q = (u 2 + v 2 ) 1 / 2 , 6 = 1 - (1 + {9 2m /9 c f)- 1 ' 2 
is the fraction of the total cluster flux that lies within 
#200 = 7*200 / dA> the angular extent of the virial radius. 

In general, all interferometers are incapable of measur- 
ing V\(0) directly, and the BIMA interferometer is con- 
strained to measure cluster visibilities at q > D/X ~ 600 
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(D is the telescope diameter). Eqn 6 shows clear that 
the shape, the size and the total integrated flux from the 
cluster affect the measured visibilities. An accurate de- 
termination of the survey mass selection function f(M) 
requires knowledge of cluster morphologies. This fact also 
means that the simple Hubble parameter h scaling that 
one derives by assuming that the total cluster SZE flux 
is the only important variable is incorrect. We discuss 
this point more in the Appendix. In future surveys it will 
be possible to use observed cluster visibilities along with 
independent mass measurements to accurately determine 
the survey selection function. In the current work we use 
hydrodynamical simulations of clusters to determine the 
survey selection. 

2.2. The BIMA Deep Survey Data 

The data are from a study of the CMB anisotropy on 
arcminute angular scales (Holzapfel ct al. 2000a). Fields 
were chosen to contain no known, bright radio point 
sources. The observations were made with the BIMA ar- 
ray with a central frequency of 28.5 GHz and ^800 MHz 
of bandwidth. The receivers are based on low-noise, high 
electron mobility transistor (HEMT) amplifiers, and are 
sensitive only to right-circularly polarized radiation. More 
details on the receivers are given in Carlstrom et al. (1996). 
Supporting observations were done with the Owens Valley 
Radio Observatory (OVRO) for aid in the detection and 
removal of radio point sources in two of the fields. 

Of the seven BIMA fields included in this analysis, three 
were observed in 1997 (the first 3 in Table 1), and the rest 
were observed in 1998. Additional observations of these 
and other CMB fields have since been obtained (Dawson 
et al. 2001, 2002). Previous observations by other teams 
(Jones et al. 1997; Richards et al. 1997) suggested temper- 
ature decrements with marginal significance in two of the 
1997 fields. However, the higher sensitivity BIMA obser- 
vations failed to confirm these detections. The four 1998 
fields are evenly distributed in right ascension, in regions 
of low IR dust emission that are free of bright optical and 
X-ray sources. Additional details of these fields are sum- 
marized in Table 1 and in Holzapfel ct al. (2000a). 

For the purposes of our analysis, only the list of visi- 
bility measurement positions (ui,Vi) and uncertainties of 
are required. We perform no further analysis of the actual 
observations; however, we describe some general details of 
the analysis from Holzapfel et al. (2000a). To remove point 
sources from the data, a long baseline (> 2.4 fcA) map was 
made for each field; five point sources were found in the 
seven fields (but see Dawson et al. 2002). Identified point 
sources were removed from the visibility data by subtract- 
ing the best fit model of the point source, modulated by 
the primary beam. 

The processed data were analyzed by a Bayesian max- 
imum likelihood method. A likelihood function was con- 
structed from the visibility correlation matrix. With this 
technique the consistency of a theory on CMB fluctuations 
could be obtained by comparison to the data. Because the 
fields are independent, the total probability is the prod- 
uct of that from individual fields. The joint likelihood for 
the combined data is expressed in terms of the flat-band 
power amplitude: Qfiat = 5Cil(l + 1)/24tt < 14.1^iK at 
95% confidence, where C; is the amplitude of the CMB 



angular power spectrum of mode I. 

Holzapfel et al. (2000a) found that, of the seven fields, 
two show excess power (Qfiat > 0); however, the confi- 
dence of a non-trivial signal is only 44% for the joint like- 
lihood analysis, and each of the fields is consistent with no 
signal at 95% confidence level. Based on this information 
together with the argument that the thermal SZE is the 
dominant source of anisotropy on this angular scale, we 
infer that no galaxy clusters were detected at or above the 
95% confidence limit. This inference is uncertain, because 
quantitative results of a sophisticated search for extended, 
negative sources (i.e. cluster SZE signatures) have not 
been reported by the authors in the three manuscripts that 
have appeared to date. Thus, we also examine the sensi- 
tivity of our limits on to the numerical significance we 
ascribe to the non-detections (see §4.3.4). 

2.3. Determining the Detection Efficiency 

The cluster mass M is the single most important factor 
in interferometric cluster SZE observations (e.g. Holder 
et al. 2000). From Eqn 5 we expect a direct correspon- 
dance between the flux from a cluster and its mass. In 
an ideal flux limited survey, this naturally leads to the 
functional form of the detection efficiency which is a step 
function: f(M > M llm ) = 1 and f(M < M Hm ) = 0, 
where Mu m is the cluster mass that corresponds to the 
minimum detectable flux of the survey. However, we ex- 
pect scatter about the total SZE flux- virial mass relation, 
because of departures from equilibrium and variations in 
cluster morphology at a fixed mass. These morphological 
variations are particularly important in an interferometric 
survey, where the observed visibility has direct sensitiv- 
ity to cluster shape and size. Therefore, we use an f(M), 
that varies continuously from f(M) = for M « M\ im 
to f(M) = 1 for M » M Hm . 

The method is that given a set of cosmological param- 
eters (here Hq, Qm and a redshift, and a position 
on the sky (a, S), we simulate observations and detections 
of an ensemble of clusters with a large range of masses. 
These mock interferometric SZE observations have char- 
acteristics such as sensitivity, u-v coverage, and primary 
beam identical to those carried out in the BIMA survey. 
We then choose the detection efficiency M\ im so that it 
reflects the mass of clusters detected at 95% confidence. 

In practice, we account for the cosmological sensitivity 
of the limiting mass through the dependence of the cluster 
flux on the angular diameter distance d,A (Eqn 5). To de- 
termine the coordinate dependence of Mu m for each BIMA 
field, we carry out sets of mock observations at 8 off-axis 
angles 9 in the range 0' < 9 < 3'. 5 and 6 azimuthal angles 
4> (although we expect azimuthal symmetry). At each of 
these field positions (9, </>), we examine the detectability of 
48 simulated clusters (observed along 3 orthogonal axes, 
§2.3.1) at redshift z = 0.5, where the cluster angular diam- 
eter distance is determined within a fiducial cosmology. To 
determine the redshift variation of the limiting mass within 
each BIMA field, we make mock, on-axis (9 = 0) observa- 
tions of these 48 clusters at different redshifts. The mass 
limit at redshift z is determined using simulated clusters 
output at that same redshift. 

For each of these mock observations, we convert the 
Compton y map determined from the hydrodynamic simu- 
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Table 1 

Basic description of the fields 





Pointing Center (J2000) 


Beam Size 


Time 


rms 


rms 


Field * 


a 


S 


("x") 


(hr) 


(MK) 


(^Jy bm _1 ) 


pl643 


16 45 11.3 


+46 24 56 


98.3x116.1 


43.1 


25.1 


191 


VLA1312 


13 12 17.4 


+42 38 05 


95.2x113.4 


35.5 


31.0 


223 


pss0030 


00 30 16.4 


+17 02 40 


99.1x115.5 


36.6 


26.6 


203 


BF0028 


00 28 04.4 


+28 23 06 


108.9x118.8 


77.6 


11.5 


99 


HDF_BIMA 


12 36 49.4 


+62 12 58 


110.9x122.0 


32.0 


18.6 


168 


BF1821 


18 21 00.0 


+59 15 00 


108.4x122.4 


43.5 


22.4 


224 


BF0658 


06 58 45.0 


+55 17 00 


105.2x148.0 


44.7 


29.3 


304 



T From Holzapfel et al. (2000a) 

*Field names are refered in Holzapfel et al. (2000a) as BDF1-7, respectively 



lations to visibilities using a Fast Fourier Transform (FFT) 
according to the prescription outlined in §2.1. Visibilities 
are then sampled from this map at the same locations in 
u-v space probed by the BIMA observations in each field. 
We add Gaussian noise to the visibility that is appropri- 
ate given the uncertainty of that same visibility in the 
BIMA observations; in fact, we scale these uncertainties 
down by a factor, which corresponds to a longer observa- 
tion with BIMA (using ay oc l/v^o&s)- This reduction in 
noise allows a more accurate estimate of the cluster mass 
corresponding to a particular detection significance (see 
Fig 1), and this scaling can be removed after the detec- 
tion. We have tested the accuracy with which one can 
introduce these longer observing times and then remove 
them after the mock detection by calculating the limiting 
mass for a wide range of observing time scale factors. Re- 
sults indicate that the limiting mass corresponding to a 
2c detection is independent of the exact scale factor used 
(assuming it is large) at the 1% level. 

2.3.1. The Hydrodynamical Simulations 

The cluster models used in the mock observations are 
from an ensemble of 48 hydrodynamical cluster simula- 
tions in four Cold Dark Matter (CDM) cosmologies: (1) 
Standard CDM (SCDM: Q M = 1, <x % = 0.6, h = 0.5, 
r = 0.5), (2) Tilted CDM (rCDM: Q M = 1, a s = 0.6, 
h = 0.5, T = 0.24), (3) Open CDM (OCDM: Q M = 0.3, 
cr 8 = 1, h = 0.8, T = 0.24), and (4) Lambda CDM 
(ACDM: n M = 0.3, n A = 0.7, a 8 = l,h= 0.8, T = 0.24). 
Here h is the Hubble parameter in units of 100 km/s/Mpc. 
The initial conditions of the simulations are Gaussian ran- 
dom fields consistent with a CDM transfer function with 
index T (Davis et al. 1985). 

The simulation scheme is particle-particle-particle- 
mesh smoothed-particle hydrodynamics (Evrard 1988, 
P3MSPH), with fixed intracluster gas mass fraction 
fiCM — 0.2 (but we correct to ficM = 0.12, consistent 
with observations for h = 0.70; Mohr et al. 1999). Ra- 
diative cooling and heat conducting are ignored. The 48 
simulated clusters have similar fractional mass resolution; 
the spatial resolution varies from 125-250 kpc. The masses 
of the final cluster sample vary by an order of magnitude. 
More detail on the simulations can be found in other anal- 
yses that have employed these same simulations (Mohr 
& Evrard 1997; Mohr et al. 1999; Mathiesen et al. 1999; 
Mathiesen & Evrard 2001). 

Following procedures described in Evrard (1990), we 



create Compton y parameter images along three orthogo- 
nal lines of sight for each cluster. For these calculations 
each cluster is imaged as though it were at a redshift 
z = 0.06. These projections serve as templates for the 
mock observations used to determine the limiting masses. 
In determining the limiting mass at a redshift z, we use the 
simulation outputs that correspond to that redshift; this 
ensures that our mass limit includes the effects of cluster 
structural evolution with redshift. 

1 4 F — i — i — th 1 1 1 1 — i — i — i — n i 1 E 




M 200 [ h 65 M g1 

Fig. 1. — The detection significance, measured as Ax 2 , of 
mock observations of the 48 simulated clusters in one BIMA field 
(BF0028). Each cluster is observed along three orthogonal lines of 
sight. The intersection of the best-fit line with the Ax 2 = 9.7 line 
(corresponding to 2a detection) determines the detection threshold 
mass M iim . Here and afterwards /i65 = /i/Q.65. 

2.3.2. Mock Cluster Detection 

We analyze the mock observations using software de- 
veloped to analyze real cluster observations (Reese et al. 
2000). For simplicity and fast computation, we model the 
ICM density distribution using a spherical (3 model (Cava- 
liere & Fusco-Femiano 1978) with (3 = 4/3. For this value 
of /3, the Fourier conjugate of the model is an exponential 
(see Eqn 6); additionally, this value is roughly consistent 
with the best fits to the SZE decrement images from the 
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simulations (< (3 >= 1.1). Furthermore, we do not in- 
clude the effects of the primary beam, and this allows us 
to bypass the FFT normally required in each fit iteration. 
This fast detection scheme should be appropriate near the 
detection threshold, where there is little information about 
the detailed shape of the cluster. 

For each cluster we fit the core radius and the central 
decrement simultaneously by minimizing the \ 2 difference 
between the fit and the observations. We do not actually 
simulate "finding" the clusters within the field; rather, we 
assume that the approximate location of each cluster can 
be obtained by examining the B1MA observation in the 
image plane. We are currently exploring other techniques, 
which would allow detection and characterization of clus- 
ters in the Fourier domain. The \ 2 difference between the 
null model and the best fit model (A% 2 = Xnull ~ A/it) 
serves as a measure of the detection significance. At a 
point on the field (9, 4>) and a redshift z, we repeat this pro- 
cedure for all 48 cluster models (viewed along three axes) 
and obtain a tight correlation between cluster mass A/200 
and detection significance A^ 2 - We determine the best- fit 
of the A% 2 — M200 relation by a downhill simplex method 
(Press et al. 1992). For 2 degrees of freedom the 95% confi- 
dence limit corresponds to A\ 2 = 9.7. The limiting mass 
Mu m (8, (f>, z) is the intersection of the best-fit line with 
the A% 2 = 9.7 line, and the detection efficiency function 
f(M) is constructed from the scatter about the best-fit 
line. Fig 1 is a plot of A% 2 versus M200 for mock obser- 
vations on field BF0028, observed on-axis (9 = <fi = 0) at 
z = 0.5. 




9(arcmin) 

Fig. 2. — The angular dependence of the limiting mass in the 
field BF0028. These limiting masses are determined using mock 
observations of simulated galaxy clusters. The mock observations 
mimic the actual observations of each BIMA field. Also shown is 
the spline, which we use to determine Mn m for arbitrary off-axis 
angle. (For the h dependence of Af ;im see the Appendix.) 

We have compared the limiting mass given by the faster 
detection algorithm described above (fixed (3 = 4/3 model 
and no primary beam) with that obtained from full /3— 
model fitting that includes the effects of the primary beam. 



The differences in the derived limiting mass Jsfl\i ra axe a/fc 
the 1% level. 

The rms scatter with respect to the A^ 2 — M corre- 
lation gives the uncertainty in the limiting mass and de- 
termines f(M). We model the detection efficiency by as- 
suming Gaussian scatter in the log A^ 2 — log M relation. 
Taking A = log A^ 2 and the rms scatter in A to be o\ , we 
can write the detection efficiency as 

f(M) = -_ / dXe n ) (7) 
\J2-ko\ Jx min 

where A TO j„ = log 9.7 (i.e. a 2a detection) and X(M) is the 
best fit log Ax 2 for mass M. 

Applying the same procedure for each deep field, we ob- 
tain Mu m (6,(f),z) and /(M) at specific values of {6,(f>,z) 
for all fields. Values of M\ im for arbitrary (6, <j>, z) re- 
quire an interpolation scheme. We take the limiting mass 
to be Mi im (9,4i,z) = M hm (6,(j),z = 0.5) x r(z), where 
Mu m (9, <fi, z = 0.5) encodes the angular dependence of the 
limiting mass at redshift z = 0.5, and r(z) is the redshift 
behavior of the mass limit at 9 = <j> = 0, normalized so 
that r(z = 0.5) = 1.0. 

— I 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 1 1 1 1 1 1 




1.5 2 2.5 3 

1+z 

Fig. 3. — The redshift dependence of the limiting mass. Points for 
each BIMA field denote the measured limiting masses at 5 redshifts 
relative to the value at z = 0.5. The best fit power law r(z) to the 
ensemble is noted in the figure. Error bars are shown only for field 
BF0028. 

To determine the angular variation of the limiting mass, 
we average the values of Mu m (d,(j),z = 0.5) over the az- 
imuthal angles <fi that are evenly spaced in the 360° circle; 
we sample off-axis angles 9 at 9 = 0^5, F0, 1'5, 2% 2'5, 
3'0 and 3'5, where the maximum 9 is close to the HWHM 
of the beam. The limiting masses at arbitrary off-axis an- 
gles are calculated use a spline. Fig 2 contains a plot of 
the angular dependence of M\ im in BF0028. 

We determine the redshift dependence r(z) using 
the limiting mass measurements at redshifts z = 
0.2, 0.5, 1.0, 1.5 and 2.3. Fig 3 contains a plot of the on- 
axis limiting mass Mn m for each field at these five differ- 
ent redshifts. A single power-law provides an adequate 
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description of the redshift behavior of the limiting mass in 
all seven fields. The scatter of the limiting masses around 
this best fit power-law is comparable to or smaller than 
the estimated uncertainties in these limiting masses. 

We have also examined the angular variation of the scat- 
ter o\ for all off-axis positions sampled (at z = 0.5). The 
variation is small (e.g. the rms variation of <j\ in the field 
BF0028 is 6.6%). There is also some redshift variation of 
a\ for each field (at 6 = (f> = 0). With a large sample of 
simulated clusters it will be possible to better determine 
the scatter in detectability <t\, but for this analysis we 
adopt a single <j\ for each field. We use this in determin- 
ing f(M) at all off-axis angles 9 and redshifts z. 

Finally, as previously noted, the cosmology dependence 
in Mu m is from the angular diameter distance to the clus- 
ters. Because we probe a wide range of cosmological mod- 
els, while our limiting mass is derived from a fiducial cos- 
mology with (Om,^a) = (0.3,0.0), we rescale Mu m to be 
consistent with the model considered. Following Haiman 
et al. (2001), the limiting mass in a given cosmology at red- 
shift z is M Um (z)h = M* lm {z)h*{d A (z)/d A \z)) 1 - 2 ^ where 
quantities with an asterisk denote that in the fiducial cos- 
mology, and cIa(z) = hd,A{z) is the Hubble-parameter- 
independent part of the angular diameter distance. We 
discuss the Hubble parameter dependence of the survey in 
the Appendix. 

3. EXPECTED YIELDS FOR ANY COSMOLOGY 

From the discussions in §2 it is clear that in calculat- 
ing the expected number of clusters (Eqn 1), Mu m , f(M) 
and the functional form of the mass function are impor- 
tant ingredients. For the mass function we adopt the fit- 
ting formula given by Jenkins et al. (2001), obtained from 
large-scale N-body simulations of structure formation: 

■jtj(M, z) = -O.SlS^a-^e-l - 61 " 1 -! 3 8 (8) 
dM ' M dM v ' 

where po is the present-day matter density and cr(M, z) 

is rms variation in overdensity on mass scale M and at 

redshift z. This fitting formula provides a more accurate 

abundance of collapsed halos than does the classical Press- 

Schechtcr mass function (Press & Schechter 1974). In both 

Jenkins and Prcss-Schechter mass functions, the form and 

amplitude of a(M) is one essential component. 

We calculate a(M) by directly integrating the power 
spectrum of density fluctuations, masked by the Fourier 
conjugate of a spherical top-hat window function W(k, R) 
over wavenumber k, where R denotes the radius of the 
spherical top-hat. We assume the initial power spectrum 
is of the form Pini(k) oc k n , which can be related to the 
power spectrum P(k, z) at any redshift after the epoch 
of matter-radiation equality by the CDM transfer func- 
tion T(k) and the growth function D{z). We use the 
transfer function given by Eisenstein & Hu (1998), and 
the baryon density parameter £Ib = 0.019 hr 2 which de- 
rives from deuterium abundance measurements (Buries & 
Tytler 1998). We write P(k,z) oc k n T 2 (k) D 2 (z) / D 2 (0) , 
where D(z) = g(fl(z))/l + z and g(Q(z)) is the growth- 
suppression factor given by Carroll et al. (1992). 

We employ a variety of strategies for normalizing o~(M). 
These include (1) matching the amplitude of the power 
spectrum on the largest scales to the constraints from cos- 
mic microwave background anisotropy measurements (e.g. 



Bunn & White 1997; Majumdar & Subrahmanyan 2000), 
(2) reproducing the local cluster abundance (e.g., Viana & 
Liddle 1999), and (3) simply adopting the normalization as 
a free parameter to be constrained using the BIMA deep 
field observations. The first approach probes the power 
spectrum on much larger scales than the last two. 

Using the Jenkins et al. (2001) expression for the mass 
function dn/dM requires adopting their halo mass defini- 
tion. They use Msoiso masses, corresponding to the mass 
enclosed within the region that has a mean overdensity of 
180 times the background density. The limiting masses we 
estimate using mock observations of hydrodynamic simu- 
lations are all M2oo's — that is, masses enclosed in regions 
with mean overdensity of 200 with respect to the critical 
density. We convert from M200 to Msoiso using the "uni- 
versal" halo density profile with a concentration parameter 
c = 5 (Navarro et al. 1997). Any differences between real 
cluster halos and our adopted model will manifest them- 
selves as a systematic error in the mass selection function 
f(M). We account for these systematics in our analysis 
(see §4.3). 

With these ingredients, we calculate the expected num- 
ber of detected clusters (N) in our BIMA fields. Rewrit- 
ing Eqn 1 to reflect our approach, the expected number of 
clusters detected in a given field is: 

m ^FeJp/J^ffMf^, , 9) 

where 9 max — 3'5, / = f(M, 9, z) is the detection effi- 
ciency and dn/dM is the cluster abundance given in Eqn 8. 
We truncate the redshift integration at some maximum 
redshift z maxi beyond which the cluster formation and evo- 
lution is uncertain. For Qm = 1-0 models, z max is set to 
be 4, and for all other models considered z max = 10. 

4. RESULTS 

We use the techniques described above to constrain cos- 
mological parameters with the BIMA deep survey fields. 
We restrict our analysis to two sets of models: (1) those 
with purely dark matter (J1m < 1 & = 0, open models 
hereafter) and (2) those with dark matter and a cosmo- 
logical constant (f^A > and w = — 1, where w = p/p 
characterizes the equation of state) that are geometrically 
flat (f^M + = 1) flat models hereafter). 

Let us denote any cosmological model to be considered 
as {f2, cr$}j. To determine the probability of consistency 
of this model with our data, we calculate the expected 
number of detected clusters (N)i for each BIMA field i. 
Suppose the detection of clusters is a Poisson random pro- 
cess, the probability of observing zero clusters in field i 
is then a s }j) = e~< N ^. The total probability of 

observing zero clusters in each of the seven fields is then 

p t otm<Ts} j ) = n i Pi- 

Below we explore the cosmological information con- 
tained in the BIMA fields themselves (§4.1), and then we 
examine constraints from the BIMA fields together with 
the local, massive cluster abundance (§4.2) and together 
with the CMB constraints on the power spectrum normal- 
ization (§4.3). We end this subsection with a discussion of 
our approach to including the uncertainties in the survey 
limiting mass, the Hubble parameter, the COBE normal- 
ization parameters, and the effects of sample variance. 
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4.1. Using Only BIMA Deep Fields 

Fig 4 contains the cosmological constraints on parame- 
ters <7g and fijv/ from the BIMA survey fields alone. Solid 
lines show the equivalent la and 3<r Gaussian probabil- 
ities for flat models, and dashed lines show similar con- 
fidence regions for the open models. We have explored 
parameter ranges 0.1 < £Im < 1-0 and 0.3 < eg < 2.0. 
These constraints are upper limits, because for fixed J7m 
it is possible to reduce the number of clusters by lower- 
ing erg. As expected, high erg can be ruled out. We fit 
the 2 a contours to give a rough constraint on erg (for 
0.1 < fl M < 1.0): <7 8 < 1.00 ft^ - 43 "" -0 - 22 (flat) and 
cr 8 < 1.01 ^-0-40^^-0.14 ( open ^ at Q5% confidence level 

In calculating these results, the uncertainties in detection 
efficiency are included (see §4.3 for a detailed discussion). 
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Fig. 4. — Constraints on Slj\,j and erg for both flat (solid line) and 
open (dashed line) universe models. The range of f!jvf considered is 
0.1 < f2jvf < 1.0, and the range of ag is 0.3 < erg < 2.0. Results 
from the Pierpaoli et al. (2001) analysis of the temperature function 
are shown (dot-dashed line) with ±ltr range. 

For comparison, dot-dashed lines in Fig 4 indicate the 
preferred values and ler confidence regions derived from 
the local temperature function of 38 massive clusters (Pier- 
paoli et al. 2001, hereafter PSW). The PSW constraints 
are only for flat cosmologies. These constraints are clearly 
much stronger than those from the BIMA survey. A larger 
solid angle and higher sensitivity SZE survey which deliv- 
ers SZE cluster detections, is required to move from upper 
limits to a preferred range of values, as in the PSW anal- 
ysis. These two surveys - a large solid angle local survey 
and a small solid angle high redshift survey - are comple- 
mentary (rather than redundant); in §4.2 we explore joint 
constraints of the deep BIMA survey and a similar but 
simpler approach to the local abundance analysis. 

We can understand the shape of the confidence region as 
follows. For fixed (Im, models of higher er s lead to larger 
(N), which are more inconsistent with the BIMA survey. 
The behavior of the probability along lines of fixed as is 
explained with the help of Fig 5. In the figure we con- 



sider three flat models and one open model with erg = 1: 
flat Qm = 0.1 (solid/red), flat VLm = 0.5 (dashed/blue), 
flat f^M = 1 (dotted/black) and open = 0.5 (dot- 
dashed/green). From top to bottom the panels contain 
the cluster redshift distribution dN/dz df2, the comoving 
differential volume element dV/dz d£l = c d\(l + z) 2 / H (z) , 
the detectable cluster number density or abundance n = 
J dM dn / dM f (M) , and the limiting detectable mass 
Mum ■ In the case of the limiting mass, we show the ratio 
of Mu m in each model to its value in the Q.m = 1 model. 
The curves in the top panel are essentially the products of 
those in the second and third panels. The behavior of the 
volume element reflects that of the angular diameter dis- 
tance, which changes only slowly beyond z ~ 1. At z > 1, 
the BIMA survey is probing an order of magnitude higher 
volume in the £Im — 0.1 model than in the Qm = 1 model. 

The behavior of the abundance evolution in each model 
is primarily a reflection of the way a(M, z — 0) is normal- 
ized: although as is set to be unity for all three models, 
the mass contained in 8/1" 1 Mpc scale (= M c i) differs for 
each model because M c i cx CIm- The result is that, at 
low mass range M > Mu m , a(M,z) of high Qm models 
is always larger. Interestingly, for masses near the detec- 
tion limit the reduction in the amplitude of a(M) for low 
Qm models at fixed as more than compensates for the 
more rapid evolution of density perturbations in the high 
J7m models. Furthermore, the steepness of the mass func- 
tion implies that contributions to the number density come 
mostly from the systems with masses near the detection 
limit. Finally, the limiting masses for low density mod- 
els are larger than that of high density models. Together, 
these effects lead to lower number densities of detectable 
clusters in the low CIm models. 

Finally, consider why the BIMA survey constrains open 
models more strongly than flat models. As also shown in 
Fig 5, the dot-dashed/green line represents an = 0.5 
open model. We see there is little difference in the volume 
element between two J7m = 0.5 models, while the differ- 
ence in abundance at high redshifts is large. This can be 
understood from the slower growth of structure in open 
models; in short, open models at fixed as produce more 
high redshift clusters than flat models and are thus more 
inconsistent with the data. 

4.2. BIMA Deep Fields Plus Local Abundance 

In this case, we use both the BIMA deep fields and a 
measurement of the local massive cluster abundance. Cal- 
culating the constraints for each Qm requires two steps: 
(1) determining the erg, which reproduces the local abun- 
dance of massive clusters above some fiducial tempera- 
ture kTx, and (2) then estimating the expected number 
of clusters (N)i for each BIMA field i. We determine as 
by requiring that n(M > M t h) matches the observed local 
cluster abundance, where M t h is the mass corresponding 
to the X-ray temperature cutoff of the sample. We adopt 
the value given by Viana & Liddle (1999), which is for clus- 
ters of X-ray temperature greater than 6.2 keV at z ~ 0.05 
(see also Henry 1997). The values of as we obtain do not 
differ much from those of Pierpaoli et al. (2001), within 
the range of 0.2 < f2jv/ < 0.7. We then calculate (N)i 
in each BIMA field and the corresponding probability for 
that cosmological model. 
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However, combining these two datasets provides very 
weak, if any, constraints on the density parameter; essen- 
tially all models considered are consistent with the data. 
Models with lower matter density are relatively more in- 
compatible, but even at the model with lowest density 
(Qm = &b), the exclusion probability is still weaker than 

1(7. 
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Fig. 5. — The behavior of the confidence regions in Fig 4 results 
from several competing effects. From top to bottom, we show the 
cluster redshift distribution, the volume element, the cluster abun- 
dance, and the limiting mass for three flat models at fixed erg = 1 
and h = 0.65: Qm = 0.1 (solid), fij^f = 0.5 (dashed), and Qm = 1 
(dotted). Also shown is an Qm = 0.5 open model (dot— dashed). 
The limiting mass panel contains the ratio of Mn m to that in the 
= 1 model. 



This method produces a lower limit on Qm because we 
are requiring the local cluster abundance to be the same 
in all cosmological models; in this scenario, those mod- 
els with slower growth of density perturbations (i.e. low 
Q.M models) produce more high redshift clusters and are 
therefore more inconsistent with our SZE survey. 



4.3. BIMA Deep Fields Plus COBE 

By combining the BIMA deep fields with a COBE nor- 
malization for the power spectrum of density fluctuations, 
it is possible to place an upper limit on VLm- The COBE 
measurements on large angular scale anisotropy can be 
used to constrain the amplitude of the power spectrum on 
the scale of the present epoch horizon (3000 ~ 4000 hr x 
Mpc). Taking the amplitude of modes on this scale to 
be 5h, we can write the power spectrum in the form 
P(k,z) = 2n 2 5 2 H {c/H Q f+ n k n T 2 {k)D 2 {z)/D 2 (0) (e.g. 
Bunn & White 1997; Eisenstein & Hu 1998). We adopt 
the fitting function for parametrization of Sh provided by 
Bunn & White (1997), assuming no contributions from pri- 
mordial gravitational waves. For the index of primordial 
power spectrum, we adopt the BOOMERANG measure- 
ment n = 0.96±£;^ (Netterfield et al. 2001). Although 
we adopt different values of Sis and Hq than those used 
by Bunn & White (1997), we do not rescale Sh^m,™) 
to account for these differences; the change in amplitude 
of 5h from this rescaling would be small compared to the 
statistical uncertainty on Sh itself. Normalizing the power 
spectrum in this way leads to as that increases with fijvf- 
Because as increases with £Im> we expect the COBE 
normalized SZE survey constraints to provide an upper 
limit on Qm- From Fig 6, we find that the 95% confidence 
upper limit on £Im in flat (open) models is Hm < 0.63 
(ft]\i < 0.73). These upper limits include uncertainties in 
the SZE flux-virial mass relation, the Hubble parameter, 
the COBE normalization and the power spectrum index, 
as we L as the effect^ of the sample variance. , — Below | we 
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The key to interpreting cluster surveys is connecting the 
observables (cluster SZE visibility in this case) to the halo 
virial mass. This generally involves a mass-observable 
scaling relation, which we term the visibility-mass relation 
- and a model for the redshift evolution of this relation. 
As has already been discussed, we use hydrodynamical 
simulations to define this visibility-mass relation, and we 
parametrize the redshift evolution of this visibility-mass 
relation by making mock observations of clusters output 
at a range of redshifts. The visibility-mass relation for an 
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intcrferometric SZE survey is similar to the better known 
mass-temperature relation in an X-ray cluster survey In 
the X-ray survey the cluster temperature is used as the 
mass predictor. 

We model these uncertainties by introducing a 40% sys- 
tematic uncertainty in the normalization of the Ax 2 -mass 
relation (see Fig 1), corresponding to a 20% uncertainty 
in the visibility-mass relation. From Eqn 5 it is clear that 
this corresponds to a systematic uncertainty in the limiting 
cluster mass (remember f(Mu m ) = 0.5) of <~ 16% (see the 
discussion in §). We model this systematic uncertainty as 
a Gaussian distribution, and marginalize over this system- 
atic by integrating the probability obtained from Eqn 9 
over the Gaussian distribution in Mu m . 

4.3.2. Hubble Parameter H and Power Spectrum 
Parameters Sh and n 

Our constraints on VIm are sensitive to the Hubble pa- 
rameter through its effect on the volume element, the am- 
plitude of a(M), the critical density and the mass limit. 
Rather than choosing a particular value of Hq, we assume 
a Gaussian distribution in Ho centered on a fiducial value 
of 65 km/s/Mpc with a lcr uncertainty of 10%. We also 
account for uncertainties in the COBE normalization am- 
plitude Sh and the primordial power spectrum index n. 
We take both the distributions of Sh and n to be Gaus- 
sian, with dispersions of 7% (Bunn & White 1997) and 
10% (Netterfield et al. 2001), respectively. The final prob- 
ability of consistency with the SZE survey is calculated by 
integrating over these Gaussian distributions. 

4.3.3. Sample Variance 

To assess the importance of sample variance on our VLm 
constraints, we follow the approach of Hu & Kravstov 
(2002). Because our survey fields are small, we approx- 
imate each field as a series of circular pillboxes, each of 
radius 3'5 and thickness Az = 0.16, arranged in increasing 
rcdshift. We calculate the covariance matrix of- between 
redshift bins i and j 

/d 3 k 
—-WiW*P(k), (10) 

where n is the number of detected clusters in the bin, b 
is the mass function weighted mean bias, D is the growth 
function and W is the Fourier transform of the pillbox 
window function at each redshift. We estimate the sample 
variance cr 2 by summing the matrix elements. We have 
only included the terms with — < 1, because other el- 
ements are negligible, as noted by Hu & Kravstov (2002). 
The exclusion probability for a given cosmology is then 

roo (M-(w>) 2 
P(n M ) oc / dfjte e -" (11) 

Jo 

up to a suitable normalization factor. It turns out, the in- 
clusion of the sample variance only weaken our constraints 
to a small extent, corresponding roughly to a 1% upward 
shift in the constraint curve in Fig 6. 

4.3.4. Field Selection and Non-detection Significance 

These fields were not chosen randomly and the selection 
criteria could potentially bias the expected number of de- 
tected clusters. Fields pl643 and VLA1312 are centered in 



directions claimed to have cluster SZE decrements (Jones 
et al. 1997; Richards et al. 1997) that have not been con- 
firmed by higher sensitivity observation (Holzapfel et al. 
2000b). The field pss0030 is in the direction of a radio- 
quiet quasar. The other 4 fields were chosen to be evenly 
distributed in right ascension, at convenient declinations, 
in regions of low dust-emission and free of bright point 
sources at 1.4 GHz, in optical and in X-ray. As pointed 
out by Dawson et al. (2001), radio point sources may trace 
large scale structure, preferentially populating clusters. If 
so, as a whole the fields may be biased to contain fewer 
than the expected number of clusters. 

Clearly, our constraints on Qm are valid only if these 
fields provide an unbiased sample of the universe. One 
could also ask whether it is surprising that no clusters 
were detected in the deepest and largest SZE survey to 
date. If the fields are unbiased, then the answer is no, if 
we live in a low matter density universe. If the fields are 
biased against finding clusters, then this would make it 
even less surprising that no clusters were detected. 

Our results also depend to some degree on our inter- 
pretation of the Holzapfel et al. (2000a) search for cluster 
SZE decrements. Specifically, we present results assum- 
ing that no clusters were detected at or above the 95% 
confidence level, whereas the authors provide no quantita- 
tive discussion. Clearly, the authors were very interested 
in finding galaxy clusters, and no firm detections were 
made (Holzapfel et al. 2000a, b). A reanalysis of the data 
searching for clusters in much the same way that we have 
searched for clusters in our mock observations would be 
ideal, and such an analysis is ongoing (Carlstrom, private 
communication) . Here we simply report the constraints on 
VLm if cluster detections in the BIMA deep fields can only 
be ruled out at the 3cr and 4cr level. In these cases, clusters 
would have to be more massive to produce these stronger 
detections, and therefore one would generally expect fewer 
detections for a given density parameter. Thus, the 95% 
confidence upper limit would weaken to Qm < 0.70 and 
Om < 0.77, respectively. In general, our results are rela- 
tively insensitive to the exact characterization of the non- 
detection significance. 

5. SUMMARY AND DISCUSSION 

We describe a systematic analysis of an interferometric 
SZE survey carried out at BIMA with the Carlstrom-Joy 
30 GHz receivers. The data were taken from observations 
of the CMB anisotropy on arcminute scale (I <~ 5500), in 
which no anisotropy and no galaxy clusters were detected 
(Holzapfel et al. 2000a, b). We use these data to study the 
allowed range of cosmological parameters as and J7m- ln 
our analysis, we consider models that are flat with a cos- 
mological constant (Om + Oa = 1) and models that are 
open with no cosmological constant. 

Our analysis hinges on our ability to calibrate the sur- 
vey sensitivity or galaxy cluster selection function. To this 
end, we use mock interferometric observations of hydro- 
dynamical cluster simulations to determine the selection 
function for each of the seven fields and as a function of 
rcdshift and position within that field. Relying on hydro- 
dynamical simulations to calibrate the survey sensitivity 
introduces significant uncertainties, which can be reduced 
once empirical calibration using large surveys that include 
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detailed followup observations are available. For this anal- 
ysis we estimate the scale of the uncertainties that come 
with using hydrodynamical cluster simulations and include 
those uncertainties in determining our constraints on the 
cosmological parameter CIm- 

With BIMA survey data alone, we obtain a 95% con- 
fidence upper limit on a 8 (0 M ): cr 8 < 1-00 fi~ / - 43 ° M " ' 22 
for flat, and cr 8 < 1.01 fi M °- 40 ° M " - 14 for open models. 
Combining the BIMA survey with external constraints on 
the power spectrum allows us to constrain Qm alone. In 
combination with the local abundance of high tempera- 
ture X-ray clusters, the BIMA data lead only to a triv- 
ially weak lower limit: Qm > &>b- In combination with 
the COBE normalization of the power spectrum, the data 
provide an upper limit on Qm- < 0.63 for flat models, 
and fijvf < 0.73 for open models at 95% confidence. In- 
cluded in the constraints on f^are 40% uncertainties in 
the cluster visibility-mass relation as well as published un- 
certainties on the Hubble parameter and power spectrum 
normalization 8r and spectral index n. We also include 
the effects of sample variance. 

These constraints are valid only if these BIMA fields 
represent an unbiased sample of the universe. The field 
selection in this survey is somewhat suspect, with four of 
the seven fields specifically chosen to exclude bright radio, 
X-ray and optical point sources. A firm interpretation of 
our results is that - if the fields are an unbiased clusters 
survey - we have shown that within the context of the 
cosmological models favored today, we expect no clusters 
to have been detected. If the survey was biased against 
finding clusters, then our constraints on D,m weaken, and 
it is even less surprising that no clusters were detected. 
Of course, if the survey were biased toward finding clus- 
ters, then finding none would imply tighter constraints on 
£Im- We eagerly await planned SZE surveys of larger and 
perhaps better chosen fractions of the universe. 

In previous work, Majumdar & Subrahmanyan (2000) 
studied constraints on Qm by comparing simulated 
SZE maps with the upper limit on the arcminute-scale 
anisotropy in CMB obtained by Subrahmanyan et al. 
(2000) using the Australia Telescope Compact Array 
(ATCA) at 8.7 GHz. Their 95% confidence upper limit 
is Q,m < 0.8 for open models. This result is clearly consis- 
tent with ours, but there are important differences in the 
analysis and the data. First, the ATCA survey has lower 
sensitivity than the BIMA survey, yielding a 95% confi- 
dence limit of Q fi a t < 25 /iK compared to the BIMA sur- 
vey constraints of Qfi a t < 14.1 /iK. Second, their focus was 
on reproducing the statistical properties of the noise in the 
sky maps, whereas we have on the non-detection of clus- 
ters. Third, our modeling technique is more sophisticated 
than theirs through its use of hydrodynamical cluster sim- 
ulations that incorporate the effects of cluster morphology 
and internal structure and the evolution of these prop- 
erties with redshift. And fourth, we model the redshift 
evolution of the cluster abundance using the results from 
the latest N-body studies of structure formation (Jenk- 
ins et al. 2001), whereas they modeled cluster abundance 



using the original Press-Schechter formalism. We believe 
that the techniques we outline in this manuscript provide 
a viable approach for analyzing the forthcoming interfer- 
ometric SZE surveys that will use much higher sensitivity 
instruments to survey larger portions of the sky. 

In short, the principle improvement of our analysis over 
previously existing analyses and estimates of survey yields 
include (1) the direct accounting for the uncertainty in 
the cluster detectability (specifically, we model the scatter 
in the detection significance-virial mass relation and we 
include the uncertainty in the normalization of the SZE 
visibility-virial mass relation), (2) the modeling of the sur- 
vey sensitivity as a function of position on the sky, (3) the 
inclusion of uncertainties on other cosmological parame- 
ters of interest, and (4) an accounting for the effects of 
sample variance on our Qm constraint. 

In addition, we have examined the h dependence of the 
survey limiting mass (see Appendix), showing that for this 
intcrferometric survey, which includes sensitivity to clus- 
ter shape and flux, the h dependence departs from the flux 
limited survey expectation Mu m oc /i~ 8 / 5 to a steeper de- 
pendence of Mu m oc h^ 11 / 6 . Similar departures would be 
expected for other interferometric surveys. 

Large solid angle, high redshift SZE cluster surveys can 
in principle provide a wealth of information on struc- 
ture evolution and cosmology. Because of its redshift- 
independent nature, the SZE is a powerful tool in studies 
of the high redshift universe. An interferometric SZE sur- 
vey has advantages such as superb angular resolution and 
the ability to separate point source contamination from 
the extended clusters (Holder et al. 2000; Kneissl et al. 
2001). Future SZE surveys will have much higher sen- 
sitivity than the BIMA survey analyzed here, and will 
be carried out with the Sunyaev-Zel'dovich effect Array 
(SZA, Holder et al. 2001), the Arcminute Microkelvin Im- 
ager (AMI, Kneissl et al. 2001), and the Arcminute Mm- 
wave interferometer for Background Anisotropy (AMiBA, 
Lo et al. 2000). These surveys are expected to detect hun- 
dreds of clusters, and these large, high redshift cluster sam- 
ples will be invaluable to cosmological studies, enabling 
tests of the mix of dark energy and dark matter in our 
universe and measurements of the nature of this dark en- 
ergy (e.g. Barbosa et al. 1996; Haiman et al. 2001; Holder 
et al. 2001; Fan & Chiueh 2001; Mohr 2001; Weller et al. 
2001; Hu & Kravstov 2002). 
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APPENDIX 

Here we examine the Hq dependence of survey sensitivity, showing how in an interferometric survey it departs sig- 
nificantly from the expected scaling for a flux limited survey. For an approximately virialized and self-similar cluster 
population, the density weighted electron temperature scales with the cluster mass: (T e ) cx {Mh) 2 ' 3 (e.g. Bryan & Nor- 
man 1998; Mohr ct al. 2000). Eqn 5 then implies that the galaxy cluster SZE flux scales as S\ cx M^h 8 / 3 . Therefore, 
the limiting mass in an SZE flux limited survey would scale as Mu m cx h~ 8 / 5 . 

In interferometric surveys, it is the cluster apparent size and morphology along with the total cluster flux that play 
important roles in determining the detectability of a cluster. The limiting mass is that mass which produces a particular 
detection significance Ax 2 , where the detection significance scales as the cluster visibility squared. From Eqn 6 it 
is clear that cluster visibility is proportional to the cluster flux, so if flux were the only important factor we would 
expect A% 2 cx Mf 00 with a = 10/3. However, the best fit Ay 2 - M 20 o relations exhibit significantly shallower slopes: 
a r~j 2.3 — 2.4 < 10/3 (at z = 0.5; see Fig 1). This underscores the sensitivity of the visibility to cluster angular size 
and shape. At a particular rcdshift, increasing cluster mass means increasing total cluster flux, but this effect is partially 
offset by the larger cluster apparent size, which reduces the signal at any particular baseline. 

There is also considerable scatter in the Ax-M 2 oo relation. We have inspected the simulated clusters that fall well below 
the best fit Ay 2 - M 20 o relation; these clusters generally appear more extended and irregular (in some cases bimodal). 
Many of these morphologically complex clusters are obviously in the midst of a major merger, and their complex shape 
causes diminution of the cluster visibility over the angular scales where the BIMA survey is most sensitive. 

To obtain the h dependence of the interferometric survey sensitivity, we need only know the slope of the Ax 2 -M 20 o 
relation and the h dependence of A\ 2 and M 2 oo- At fixed ICM temperature (T e ), M 2 oo cx h^ 1 and A\ 2 cx h 2 . Including 
the h dependence in the Ax 2 -M 20 o relation leads to (A\ 2 h~ 2 ) cx (M 2 ooh) a . Keeping the required detection significance 
Ax 2 fixed provides h scaling: M/ im cx /i _ ( 2+a )/ a . For a = 10/3 which is derived from assuming only flux is important 
in determining the cluster detection significance, we arrive at Mu m cx /i~ 8//5 , consistent with expectations for a flux 
limited survey. Taking the typical measured slope a — 2.4, we calculate the true h sensitivity of the limiting mass to be 
Mi im cx hT ^ . The shallower the Ax 2 -M 20 o relation, the larger shift in Mu m required to offset a change in h. 



